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Abstract —This paper presents the performance analysis of the 
self-interference cancellation (SIC) methods in full-duplex large- 
scale multiple-input multiple-output systems. To mitigate self¬ 
interference (SI), we assume that the full duplex-base station 
(BS) uses Sl-subtraction or spatial suppression. Analytical and 
numerical results confirm that the Sl-subtraction outperforms 
the spatial suppression for SIC in a perfect channel estimation 
case. It is also concluded that the uplink and overall ergodic rates 
performance of the spatial suppression is respectively better than 
those of the Sl-subtraction in a imperfect channel estimation case 
under a given system constraint such as uplink/downlink sum 
rates and the total transmit power at the BS. 

Index Terms —Full duplex radios, large-scale MIMO, massive 
MIMO, self-interference cancellation. 

1. Introduction 

By the year 2020, it is expected that mobile data traffic 
will have grown a thousand-fold (lOOOx). To support such 
tremendous data traffic, researchers should consider in a next 
generation communication system such key issues to increas¬ 
ing network capacity (bps/area) as high spectral efficiency 
(bps/Hz), spectrum extension (Hz/cell), and high network 
density (cell/area) should be jointly considered. Over the last 
few decades, multiple-input multiple-output (MIMO) wireless 
communication techniques have evolved for high spectral 
efficiency 1^. Eull-duplex transmission has also been studied 
for double spectrum extension 0, 0 Eor high network 
density, researchers have considered deploying a large number 
of small cells. 

Recently, to maximize spectral efficiency and to conserve 
energy, the researchers in 0-0 have proposed large-scale 
MIMO systems. They considered simple linear precoders 
and filters to mitigate the interference that arise from large- 
scale MIMO transmission. The authors in 0,0 studied the 
downlink performance of maximum ratio transmission (MRT) 
and zero-forcing (ZE) precoders while assuming perfect and 
imperfect channel estimations at the base station (BS). The au¬ 
thors in |[7| investigated the uplink performance of maximum 
ratio combining (MRC), ZE and minimum mean square error 
(MMSE) filters while assuming perfect and imperfect channel 
estimations at the BS. They showed that the transmit energy of 
the users could be conserved by the power-scaling law when 
the BS with a large-scale array serves a small number of users 


that equipped with one antenna. The authors in analyzed 
both downlink and uplink performances of the cell-boundary 
users (including the performance analysis of a large number of 
users) considering the impact of the low transmit or received 
power. They also proposed downlink precoding normalization 
methods and transceiver mode selection algorithms when the 
BS uses the MRT/ZE precoder at downlink and the MRC/ZE 
filter at uplink. 

Another relevant technology is the full duplex transmission 
that can double the spectrum bandwidth compared to half¬ 
duplex transmission, such as time division duplex (TDD) and 
frequency division duplex (EDD). In contrast to the half duplex 
transmission, the self-interference (SI) which is a transmitted 
downlink signal and directly received at the BS while receiving 
a uplink signal should be mitigated. The authors in 0 com¬ 
pared the self-interference cancellation (SIC) by Sl-subtraction 
and that by spatial suppression in a MIMO-relay system with a 
small number of antennas (up to 4 antennas) at the BS and the 
user. They showed the performance of the spatial suppression 
was better than that of the Sl-subtraction but depended on the 
rank of the SI channel and that it required additional antennas. 
In a practical system with a small number of antennas at 
the BS, the Sl-subtraction is a widely considered the SIC 
technique due to the limitation of the requirement of additional 
antennas. The researchers in pQ| prototyped SI-subtraction- 
based algorithm on a software-defined radio (SDR) platform 
in real time. They combined a dual-polarization antenna-based 
analog cancellation with the Sl-subtraction-based algorithm 
and achieved 1.9 times higher spectral efficiency than the half 
duplex transmission. 

There are some integrated techniques for high spectral effi¬ 
ciency, spectrum extension, and high network density such as a 
large-scale MIMO system with a small cell technique and full 
duplex transmission with a large-scale antenna array, which 
can support huge data traffic demands. The researchers in GD 
presented a real time three-dimensional hybrid beamforming 
for next generation communication systems, which is equipped 
with the large-scale antenna array at a small cell. (D) 
proposed a compressed channel feedback for a highly cor¬ 
related channel that is yielded from compact antenna arrays 
at the small BS. Similarly, CD proposed a novel channel 
feedback for a radio frequency lens-embedded large-scale 
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Fig. 1: A block diagram of a full duplex multi-user large-scale 
MIMO system. 


MIMO system, and analyzed the performance at a small cell. 
The authors in (ID analyzed the performance of the nonlinear 
SIC techniques for a full duplex system at a small cell. |T^ 
proposed the spatial suppression-based algorithms such as 
extended ZF and extended MMSE precoders for full duplex 
large-scale MIMO systems. Since their approximations of the 
downlink/uplink sum rate are based on an interference-free 
channel, their analysis lack certitude. Moreover, the authors 
ignored the impact of imperfect channel estimation-the key 
issue in SIC problem. 

In this paper, we present a performance analysis of the SI- 
subtraction and the spatial suppression in full duplex large- 
scale MIMO systems. The spatial suppression could be a good 
option of the SIC in large-scale MIMO systems because many 
additional antennas can be utilized. From the analysis and the 
numerical results, we investigate the SIC algorithm which is 
better. 


II. System Model 

In this section, we introduce the basic notation used in this 
paper and the full duplex large-scale MIMO system^ 

A. System Model: Full Duplex BS with Large-Scale Arrays 

Consider a full duplex multi-user MIMO channel, as illus¬ 
trated in Fig. with M(> N K) transmit antennas and 
N receive antennas at the BS (a total number of antennas 
at the BS is L, i.e., L = M N) that serves K downlink 
users and K uplink users equipped with one antenna for 
each. This asymmetric antenna system ensures the small and 
low cost BS by reducing radio frequency chains that include 
power amplifiers (N > M K) Sit a. transmitter or low-noise 
amplifiers (M > -1- AT) at a receiver in the fixed size of total 

number of radio frequency chains (i.e., L). We assume the BS 
has more transmit antennas than receive antennas to support 
the asymmetric data traffic load which the downlink has more 

^Throughout this paper, we use upper and lower case boldfaces to describe 
matrix A and vector a, respectively. We denote the inverse, transpose and the 
Hermitian of matrix A by A~^, A^, and A*, respectively. 


offered than the uplink in practice. We also assume that there 
is no user-to-user interference where the path-losses between 
the users are sufficiently large. The free-space propagation gain 
^si is considered at the SI channel. The total transmit signal- 
to-noise ratio (SNR) at the BS, the received downlink SNR at 
the users, the received uplink SNR at the BS, and the received 
SI SNR at the BS are respectively denoted by pt{= Pt/Pn), 
/3dl(= Ptl3\jE), PuL’ and p^,{= ptPsi), where Pt, Pn, and /3 ue 
are the total transmit power, the noise power of the additive 
white Gaussian noise of the channel, and the propagation gain 
between the BS and the user respectively]^ We assume that 
the SI, downlink and uplink channels are Rayleigh channel 
modeled by the x M matrix the AT x M matrix F^dl, 
and the N X K matrix respectively, the elements of 

which are i.i.d. (independent identically distributed) complex 
Gaussian random variables with zero mean and unit variance. 

Let Qk denote the k-th column vector of transmit precoding 
matrix G and Sk represent the transmit symbol for the k-th 
user at downlink. Similarly, let Wk denote the k-th column 
vector of receive combining matrix W and Xk represent the 
transmit symbol for the k-th user at uplink. Also, let Uk be 
the additive white Gaussian noise vector of the k-th user. The 
received signal at the k-th user is then expressed by 

K 

Vk = v^^DL,fe5feSfe + k9ese +n^uk (1) 

---^ i=iy^k 

desired signal v ^ > 

user-interference 

where Adl,/c denotes the the k-th column vector of Also, 
the received signal for the k-th user at the BS is expressed by 

K 

rk =/p^wlh\]i^,kXk + L! '^k^vh/Xi 

desired'signal ^ 

user-interference 

+ \/Jj^cwlHsiGs+w'ln^^k ( 2 ) 

'-V-" 

self-interference 


where Aul,/c and ctanc denotes the k-th column vector of JTul 
and the level of the passive analog SIC. To satisfy a power 
constraint, the normalized transmit beamforming vectors (the 
columns of the precoding matrix) with vector normalization 
are given as gt = fk/{VK\\fk\\) l[^- 
From 0 and 0, we can derive the SINR of the k-th user 
at downlink and uplink, respectively, as given by 
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^There is a trade-off between a downlink rate and a uplink rate according 
to Pt if the values of the propagation gains are fixed. To simplify this effect 
in the expressions of the signal-to-interference-noise-ratio (SINR), we define 
the total transmit, received downlink, uplink, and SI SNRs. We also assume 
the distances from the BS to all users are the same and there are no shadow 
fading for a simplicity. 
































where O = \ \w^HsiG\\‘^. 

B. Precoding/Receive Combining Matrix Design and Self- 
Interference Cancellation Methods 

Let ^dl, H\]\^ and ^si be the estimated channel matrices 
of and Since the full duplex-BS always 

perfectly knows its precoding/receive matrix and transmitted 
downlink symbols, it can mitigate SI by directly subtracting 
of which the SIC is referred to as the SI- 
subtraction. The received signal after Sl-subtraction is 

J’stt.fe = rk- \/pJc^wlHsiGs. 

A ZF precoder for downlink is assumed at the BS, in the cases 
of the the full duplex-BS with the SI-subtraction and without 
the SIC: 

Fzf = H*r,i^{HuLH*uLr^. 

From the antenna configuration of M > N-\-K, the additional 
antennas for the spatial suppression can be utilized. The spatial 
suppression-based algorithms for a large-scale MIMO system 
was proposed in IK. The spatial suppression-based precoder 
is given by 

rp _ IT* / TT TT^ 1 TT _ DL 

^ sps — ^extV^ext^extJ ^ ^ ext — £y 

To eliminate user-interference signals at uplink, we use fol¬ 
lowing ZF receive combining matrix in both transceiver with 
the Sl-subtraction/spatial suppression and without the SIC: 

W = 

C Estimated Channel Model 

We assume pilot symbols are orthogonally transmitted under 
TDD-based transmission. The number of pilot symbols is K 
and N at the downlink/uplink channel and the SI channel. 
Since estimation errors are independent with original channels, 
estimated channels can be modeled by estimated error matrices 
{E) that have i.i.d. CA/’(0,e^) elements as follows: 

^DL/UL/SI = ^DL/UL/SI + -^dl/ul/si* (^) 

We assume that MMSE-based channel estimation is used at 
the BS Q. 

1) SI Channel Estimation Error: If the hardware of the full 
duplex system is perfect, the MMSE-based estimation error 
of the SI channel could be close to zero due to high received 
pilot power. In practice, the estimation error of the SI channel 
yields residual SI that impact on the significant performance 
degradation of the full duplex systems which error is imperfect 
because of hardware impairments such as the nonlinearity of 
the power amplifier, a quantization error, and phase noise (TT) 
Thus, it is difficult to model the estimated SI channel while 
considering the hardware impairments. The variance of the 
estimation error is easily modeled as the fixed value of the 
normalized-MSE (NMSE) between the estimated channel and 
the original channel, which depends on the performance of the 


estimation algorithm and the specification of the hardware. The 
variance of the estimation error of the SI channel is given by 

= NMSE. 

2) Downlink/uplink Channel Estimation Error: In the case 
of the downlink/uplink channel estimation, we assume there 
are no error terms from the hardware impairments since the 
distance between the BS and the users are sufficiently large, 
i.e., the path-loss is large, which means the errors from 
the hardware impairments become relatively smaller than the 
additive white Gaussian noise. Thus, the MMSE estimation 
error can be considered Q. The variance of the estimation 
error for the downlink/uplink channel is given by 

2 ^ ^UE ... 

KpJ^ + 1 

where is the transmit SNR at the user. Note that any 
estimation algorithm dependent on the transmit pilot power 
can be modeled similarly. 

III. Performance Analysis 

In this section we present the performance analysis of the 
SI-subtraction and the spatial suppression in full duplex large- 
scale MIMO systems. 


A. Perfect Channel Estimation 

Let = E{0} be expected SI power. The Sl-subtraction 
and the spatial suppression transceivers with perfect channel 
estimation promise = 0. If there is no SIC at the BS (a 

a nc ~ 

0), = E { I jtu/e I p } because Hsi and G are independent | [f^ . 

Thanks to the approximation methods of j^, we can derive 
the approximations of the ergodic sum rate of the full duplex 
transceiver without the SIC, with the Sl-subtraction, and with 
the spatial suppression, respectively, as follows: 
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where the approximations of the ergodic sum rate at downlink 
result from E | m. ^ ,12 \ = M — K and E | m. ^ ,12 \ = 
M — {N A K) A I as well as those at uplink are derived from 
E 1 1|^^||2 1 = N — K Al. The approximations of the ergodic 
sum rate of the half duplex with ZE precoding and ZE receive 
combining is the same as ^T^stt l[^- 













B. Imperfect Channel Estimation 


In this section, we analyze uplink SINR, which could 
be a more meaningful element than downlink SINR in the 
performance perspective at a full duplex system because of 
SI. From Q and the SINR of the k-th user at uplink in a 
full duplex system with imperfect channel estimation and the 
ZF receiver at the BS is given by 


7zf,/c ~ 




UL,fe| 


PUL^ + Psi^+ IKIP 
Pm. 


PmEt, |<euL,.|2 + ^fi+|K||2 
where euL,/c denotes the k-th column vector of E 


UL-. 


( 10 ) 


/ = 


+ KWP = 1, 

and = 0, respectively. 

1) SINR of a full duplex receiver without the SIC at 
uplink: By using the approximation methdos of 0, 0, the 
approximation of the uplink SINR without the SIC is derived 
from (10) and Cl = E [\\wk\\‘^} as follows: 

„ K,l{N-K) 

■ 2Kp.^ + ^{Kp^ + l) + l' 


2) SINR of SI-Subtraction at Uplink: The expected SI 
power of the full duplex transceiver with the SI-subtraction 
is 


at. = E{|K(/fsi-^si)G|P} 

= ^{\\wlEGf] 

^=72e{|K||2} (12) 


where (a) results from |T^ since E and G are independent. 
From ( p^ and {V2), we can derive the approximation of the 
uplink SINR with the Sl-subtraction, which is referred as 
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7stt,/c 
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(13) 


3) SINR of Spatial Suppression at Uplink: Since H^i and 
Gsps are not independent in the case of the spatial suppression, 
we should consider the best and the worst estimation error to 
make the independent condition. In the case of the best channel 
estimation error (e^j ^ 0), the estimated SI channel tends to 
the original SI channel {Hsi Hsi) so E and Gsps tend to be 
independent. Thus, the expected SI power of the full duplex 
transceiver with the spatial suppression (the best case) is given 
by 

(isps(e,j^0) =E{||t«fei?siGsps|P} 

= e{|K(Hsi-£;)Gsps|P} 

= E{\\wlEGs,s\f} 

= e2E{|K||2} (14) 

where .ffsi^sps = 0. Similarly, in the case of the worst channel 
estimation error (e^^ ^ oo), the estimated SI channel tends 
to the channel estimation error matrix (Hsi E) and .ffext 


tends to 


^DL 

E 


so Hs\ and Gsps tend to be independent. Thus, 

the expected SI power of the full duplex transceiver with the 
spatial suppression (the worst case) is given by 

iisps(e,. ^oo) = E { I \wlHsiG,^s 11" } = E {I 112 } . (15) 

From and the approximation of the expected SI 
power of the full duplex transceiver with the spatial suppres¬ 
sion is expressed by 

.E{|K|| 2 } 


a 


sps 


1/el + I 


This approximation approach is a simple application of the 
harmonic mean. Then, the approximation of the uplink SINR 
with the spatial suppression are given by 

KpiiN-K) 

Tsns k — 7 \ • (lb) 

“^^Pm + ^ (l/4+l) i^PvL + 1) + 1 

C. Comparison of SI-Subtraction and Spatial Suppression 

From ^ and ([^, the SI-subtraction outperforms the spatial 
suppression in perfect channel estimation. Next, we compare 
the expected SI power of the spatial suppression and the SI- 
subtraction as follows: 
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= e2E{|K||2} =att 

where the inequality is obtained by using Cauchy-Schwarz’ 
inequality. C denotes EEi \gk,m?\El^Ei^j\ and Cde¬ 

notes EEi Yf=i E| 5 fe,m|^E|£’*,„£’ij| where Wij, gij, and 
Eij are the (i, j)-th element of W, Gsps, and E, respectively. 
From we conclude that the ergodic performance of the 
spatial suppression is better than those of the Sl-subtraction 
in imperfect channel estimation. 


IV. Numerical Results 

For numerical comparisons, we assume that the full duplex- 
BS has a total of 84 antennas (M = 64, = 20) and serves 10 

users for each downlink and uplink. Simulation parameters are 
set as Pul = 10 /^si = —40 dB, and /3 ue = —80 dB, which 
are reproduced from the typical values of the path-loss at the 
reference meter (1 m) and those at small cell environments. 
We also assume that (Tacn = 40 dB and NMSE = 0.2. 


A. Proposed Analysis and Simulation Results 

We first compare the proposed analysis and the simulation 
results. Figure shows that the results from (|7]), ([^, and 0 
are approximately the same as the ergodic downlink/uplink 
sum rate of the full duplex-BS without the SIC and with 
the SI-subtraction/spatial suppression where perfect channel 














Fig. 2: downlink/uplink sum rate vs. received downlink SNR 
with perfect channel estimation. 



Fig. 3: uplink sum rate vs. received SI SNR with imperfect 
channel estimation 


estimation is assumed at the BS. Since the SI is perfectly 
canceled out due to perfect knowledge of the SI channel, 
both uplink performance of the spatial suppression and the 
SI- subtraction are the same which use the same ZF receiver, 
i.e., there are a performance gap only in the downlink case. 
Additional spatial domain (A^) is needed to cancel out the 
SI channel at the full duplex-BS with the spatial suppression, 
so we confirm that the SI-subtraction outperforms the spatial 
suppression in perfect channel estimation. 

Similarly, Fig. shows that the approximation forms, ([n), 
and ( p^ , close to the ergodic uplink sum rate of the full 
duplex-BS without the SIC and with the SI-subtraction/spatial 
suppression where imperfect channel estimation is assumed. 
The legend. Passive ANC, indicates the ergodic sum rate of the 
full duplex-BS only with a passive analog SIC and without any 
digital SIC. Note that these approximations are more tighter 
at the high received SI SNR regime, which are based on 
the approximation method for the ergodic sum rate in large- 
scale MIMO systems at the high SNR regime in |[^. From 
the numerical results and the analysis, we can conclude that 
the spatial suppression technique could mitigate both SI term 
and estimation error term due to the correlation between the 
spatial suppression precoder and the estimation error of the SI 
channel. 

B. Simulation Results for a Correlated Channel 

Next, we compare the Sl-subtraction and the spatial sup¬ 
pression in a highly correlated channel. To support huge 
data traffic demands, the full duplex BS is assumed to be 
located in a small cell, giving it compact antenna arrays with 
a high channel correlation|^ For numerical comparisons, the 
SI channel is a Rician channel modeled by the x M 
matrix with the line-of-sight (LOS) channel variance cr|j 

^ A SDR platform for a small cell with compact multiple layers of a uniform 
linear array was presented in 0- 


and Rician factor k. Let iJxx and be spatial correlation 
matrices of the transmit antennas and the receive antennas at 
the BS. From Jakes’ model, the element of i?Tx and Brx 
is = Jo (27rJ™^/A) where is the distance 

between the i-th array and the j-th array of the transmit/receive 
antenna at the BS, and Jo and A are a zero-order Bessel 
function of the first kind and a carrier wavelength, respectively. 
Spatially-correlated DL, UL, and SI MIMO channels can be 
respectively modeled as: 


rr _ rr TJ _ 

il DL — ^ DLiid-^TX 5 UL — -«RX ^ ULno ^ 


^si = R. 
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+ 1 


C^Sl 




SllID 
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1/2 

TX 


where Hulud and JTsinD are respectively i.i.d. Rayleigh 

DL, UL, and SI MIMO channels. The system operates at 
2.1 GHz of the carrier frequency and imperfect channel 
estimation is assumed at the BS. 

In Fig. 1^ we consider a highly correlated SI channel with 
LOS components. Thus, a Rician channel is assumed with 
K = 1 and (Jsi = 1, and each path-loss between the i-th receive 
array and the j-th transmit array is modeled by free-space 
path-loss according to dij and the minimum value of dij is | 
with a uniform linear array. The result shows that the uplink 
sum rate of the spatial suppression outperforms that of the SI- 
subtraction while the downlink sum rate of the Sl-subtraction 
better than that of the spatial suppression. 

The more the downlink sum rates are increased, the more 
the uplink sum rates are degraded because high transmit 
power produce high SI power in a full duplex system. We 
should consider the performances at downlink and uplink, 
simultaneously. Figure]^ shows that the overall sum rate of the 
spatial suppression outperforms that of the Sl-subtraction for 
the given constraint of the uplink sum rate (e.g., 10 bps/Hz). 
Furthermore, the overall sum rate of the spatial suppression 



























































Fig. 4: downlink/uplink sum rate vs. received downlink SNR 
with imperfect channel estimation where channels are corre¬ 
lated. More realistic parameters are used in simulation. The 
legends, UL and DL, indicate uplink and downlink, respec¬ 
tively. 


better than that of the SI-subtraction for the given constraint 
for the downlink sum rate where the uplink sum rates are 
guaranteed (e.g., at low pdl (< 15 dB)). We summarize these 
conclusions in Table U 

V. Conclusion 

This paper has presented the performance analysis of the 
Sl-subtraction and the spatial suppression in full duplex large- 
scale MIMO systems. The analysis and the numerical results 
have revealed the SI-subtraction outperforms the spatial sup¬ 
pression with perfect channel estimation because the spatial 
suppression requires additional antennas for null projection. 
We have also concluded that the ergodic performance of the 
spatial suppression is better than those of the Sl-subtraction at 
uplink with imperfect channel estimation due to the correlation 
between the spatial suppression precoder and the estimation 
error of the SI channel. We provided an insight into which 
SIC algorithm is better under a given system constraint such 
as the downlink/uplink sum rate and the total transmit power 
at the BS. In future work, we will investigate the impact of 
the hardware impairments, user-to-user interference, and pilot 
overhead issues. 
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TABLE I: Desired SIC techniques for a given system con¬ 
straint 


Given System Constraint 

Desired SIC Technique 

Downlink Sum Rate 

Spatial Suppression 

Uplink Sum Rate 

Spatial Suppression 

Total Transmit Power 

Sl-Subtraction 
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